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ABSTRACT 

Aims. Some theoretical scenarios suggest the formation of brown dwarfs as ejected stellar embryos in star-forming clusters. Such 
a formation mechanism can result in different spatial distributions of stars and substellar objects. We aim to investigate the spatial 
structure of stellar and substellar objects in two well sampled and nearby embedded clusters, namely IC348 and the Orion Trapezium 
Cluster (OTC) to test this hypothesis. 

Methods. Deep near-infrared K-band data complete enough to sample the substellar population in IC348 and OTC are obtained from 
the literature. The spatial distribution of the K-band point sources is analysed using the Minimum Spanning Tree (MST) method. The 
Q parameter and the spanning trees are evaluated for stellar and substellar objects as a function of cluster core radius R f . 
Results. The stellar population in both IC348 and OTC display a clustered distribution whereas the substellar population is distributed 
homogeneously in space within twice the cluster core radius. Although the substellar objects do not appear to be bound by the cluster 
potential well, they are still within the limits of the cluster and not significantly displaced from their birth sites. 
Conclusions. The spatially homogeneous distribution of substellar objects is best explained by assuming higher initial velocities, 
distributed in a random manner and going through multiple interactions. The overall spatial coincidence of these objects with the 
cluster locations can be understood if these objects are nevertheless travelling slowly enough so as to feel the gravitational influence 
of the cluster. The observations support the scenario of formation of substellar objects as "ejected stellar embryos". Higher ejection 
velocities are necessary but net spatial displacements may not be necessary to explain the observational data. 

Key words. Stars:low-mass, brown dwarfs; Stars: formation ; ISM: kinematics and dynamics; Methods: statistical 



1. Introduction 

Observations of embedded clusters display structure and mass 
segregation generally thought to reflect the structure of molec- 
ular clouds from which they were born an d the gravitational 
potential of the cluster (lLada & L ada 2003). Massive stars are 
found to be located at the centers of the clusters and the ra- 
dial distribution of cluste r members roughly follow King models 
or power law variations (Hillenbr and & Har tmann 1998). These 
findings are generally based on observations sampling the stel- 
lar content of the clusters. With the advent of sensitive infrared 
detectors, the substellar content has been effectively unveiled 
allowing studies of brown dwarfs and planetary mass objects. 
In the rest of this paper we shall refer to Brown Dwarfs and 
free-floating planetary mass objects as BDs without distinction. 
The BDs are known to be particularly luminous in their early 
stages owing to which hundreds of such ob jects are detected in 
near-infrared surv eys of embedded clusters dMuench etaiT2 003: 
iLucas et alj|2005l) . 

The formation of substellar objects can be explained throug h 
at least five alternative mechanisms (W hitworth et al.l [2007). 
One of the important mechanisms, the so-called "ejected em- 
bryos" scenario, advocates formation of BDs as a result of 
premature ejection of protostellar embryos from multiple sys- 
tems. This "ejection scenario" attributes higher velocity dis- 
persion an d spatial spread to BD s in comparision to stel - 
lar objects dReipurth & Clarkg|2001l: iKroupa & Bouvierf l2003). 
Observations of the radial velocities of BDs and stars in Taurus 
and Chameleon star forming regions indicate that there is no 



significant differen ce in the velocity d ispersion between stars 
and brown dwarfs dLuhman et al .1120071) . Young stars and BDs 
are also thought to be homogeneously mixed in the Trapezium 
cluster and Taurus molecular cloud. Much of these results are 
obtained by studying a few objects in localized regions of the 
cluster or in a region like Taurus where star formation occurs 
in relat ive isolation. In a dee p near-infrared study of the IC348 
cluster, Muenc het al.l (120031) analysed the radial distribution of 
point sources which indicated a relatively homogeneous distri- 
bution of BDs in space compared to the stellar population that is 
bound to the cluster potential. 

The structure analysis of embedded clusters has gained im- 
portance in recent years, since numerical simulations of the star 
formation process on the scales of molecular clouds are avail- 
able (e.g. iBateet all 120031) . which can be tested against obser- 
vational results. Further, to facilitate the structure analysis of 
observed data, quantitative methods of statistical analysis have 
also been developed. The nearest neigh bour method and the min- 
imum spanning tree (MST) method dCartw right & Whitworth 
2004) are two statistical methods that can be applied on obser- 
vational data to analyse and quantify the physical parameters 
of clusters based on structure analysis. The structures of some 
star forming clusters have been analysed using the MST m ethod 
dCartwright & Whitworthll2004tlSchmeia & Klessenll2006l) . 

IC348 and the Orion Trapezium cluster (hereafter OTC) are 
two nearby well-studied star forming regions for which deep 
near-infrared data are available. The OTC in particular has been 
the target of several deep surveys that aimed to unveil the sam- 
ple of BDs down to planetary mass objects. In both, IC348 and 
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OTC, the BD population is well sampled. These are also richly 
populated clusters where a statistically significant sample of sub- 
stellar objects are catalogued. In this letter we will analyse the 
structure of both IC348 and OTC to study the variations for the 
stellar and BD populations. 



2. Data sets and analysis methods 

Near-infr ared point source cata logs for the OTC were obtained 
from the IMuench et"ai] d2002l) and iLucas et all d2005l) st udies 
via the on-line VizieR database. The IMuench et al.l ([2002) data 
cover a region centered on the Trapezium cluster to a 10cr de- 
te ction limit of K - 18.1 mag and inclu de all the sources f ound 
bv lHillenbrand & Carpenter! feoOOt) . The lLucas etail (120051) data 
obtained through FLAMINGOS on GEMINI south telescope 
cover a region adjacent to the boundary of the Trapezium cluster 
which suffers from relatively little contamination from the K- 
band nebulosity. We combined both the catalogs taking care of 
overlaps that resulted in a coverage of ~ 0° 1 x 0° 1 . In the over- 
lapping regions between the two surv eys we chose to use the 
magnitudes from IMuench et all d2002l) . Data for the IC348 re- 
gion were obtained from the FLAMINGOS survey described by 
IMuench et al.l d2003l) . The data cover an ~21'x21' area centered 
on the IC348 cluster to a depth of K~17 mag. 

Based on comparisons to pre-main sequence evolution- 
ary tra cks (|Lucas et at] 20051) and luminosity function mod- 
elling (IMuench et al.l 120021) . in the OTC sources brighter than 
K=14 mag can be considered stars and sources with K>14 can 
be considered BDs. Similarly, in IC348, sources fainter than 
K=15 mag should be complete to O.O4M brown dwarfs seen 
through an extinction of A,, ~7. These estimations take into ac- 
count the effects of extinction and nebulosity, hence they repre- 
sent a sample not contaminated by such effects. We will use this 
classification in the following analysis to distinguish stars and 
BDs. 

Control field data were obtained from the 2MASS point 
source catalog. The control fields were selected by looking at 
large area (5°) IRAS 100 /mi maps and choosing a region of very 
low emission. Point sources within a circle of radius 15' were 
extracted with the criteria of photometric quality flags better than 
"D" and without confusion (ccflag=0). The Orion control field 
is centered at 05 /, 24 m -07°20'. In the case of IC348, the off- 
cloud control field is located at 03 /! 36 m +35°00'; another field 
on the Perseus molecular cloud centered at 03 /! 39 m +31°20' was 
also examined. Although the 2MASS survey 100% complete- 
ness limit is at K=14 mag, several reliable detections of K=15- 
17 mag are also found in the point source catalog at a lower <x 
level. 



We use the positions and magnitudes to analyse the spatial 
distribution of the objects by means of a minimum spanning tree, 
the unique set of straight lines ("edges") connecting a given set 
of points without closed l oops, such that the s um of the edge 
lengths is a minimum (e.g lGower & Rossl|l96 9). We derive the 
mean edge length ^mst as well as the parameter Q - £mst/s by 
combining the normalised mean edge l ength and the normalised 
correla tion length s as described in [Cartwright & Whitworth 
d2004 and lSchmeia & Klessenl d2006l) . The Q parameter quan- 
tifies the cluster structure and distinguishes between a centrally 
condensated cluster (Q > 0.8) and a cluster showing fractal sub- 
structure (Q < 0.8), while ^mst can be seen as the mean distance 
between two neighbou ring o bjects. We construct the MST using 
the algorithm of lPriml d!957l) implemented in an IDL routine. 



3. Results 

The MST method was applied on the photometric data to analyse 
the following: 

a) The variation of the Q parameter as a function of source 
magnitude (therefore source mass) selected to construct the tree. 
At every magnitude (intervals of 0.2 mag), Q was evaluated by 
considering all sources brighter than the given magnitude. This 
variation was studied for the cluster field and the control fields. 

b) The structure of the minimum spanning tree for stars and 
BDs separately. The resultant mean edge lengths ^mst along 
with their standard deviation and minimum and maximum val- 
ues were obtained for each tree. Projected separations less than 
150 AU (indicating source binarity) are neglected in computing 
^mst to remove their effect on the overall distribution. Further, 
the edge lenghts were examined inside and outside a radius that 
define the cluster core boundaries. 

In Figures 1 and 2 we show these results for IC348 and OTC. 
In both figures, panel a) shows the variation of Q with magni- 
tude, while panels b) and c) show the tree structure for stellar 
and substellar sources respectively. The dotted circles represent 
boundary limits of the cluster as defined by radial profile mod- 
elling. The ^mst values for each zone enclosed by the bound- 
aries, along with the statistics are listed in Table 1 . 

It can be seen from Figs, la and 2a that the Q values peak 
at brighter magnitudes and show a turn-over at K=14 mag for 
IC348 and K=ll mag for OTC decreasing smoothly towards 
fainter magnitudes. In comparision, the data from the control 
fields do not show significant variations of the Q paramater with 
magnitude. They yield values in the range of Q « 0.72, the 
value expected for a random distribution. The peak Q values 
for IC348 and OTC are 0.9 and 0.94, respectively, indicating 
centrally condensed configura tions. Following the findings of 
ICartwright & Whitworthl d2004l) (condensation for Q > 0.8, frac- 
tral substructure for Q < 0.8), these results imply that in IC348 
sources fainter than K=15.5 mag (the point where Q touches 
0.8 in Fig. la) are distributed in space with a fractal structure 
while brigther sources are more centrally condensed. The same 
intrepretation in OTC suggests a general central condensation 
that is only weak for the fainter sources. 

The structure of the MST in Figs, lb and 2b display clusters 
at the center of the field while Figs, lc and 2c display a ho- 
mogeneous distribution of sources. The clustering can be seen 
as branches of smaller edge lengths at the center spreading into 
branches of longer edge lengths away from the center. It can 
also be noted that the clusters in Figs, lb and 2b coincide well 
with the B star in IC348 and with the Trapezium stars in the 
OTC cluster, shown with filled star symbols. The fainter sources 
do not show detectable variations in the tree structure across 
the field. These results are demonstrated more quantitatively in 
Table 1 . The linear lengths in Table 1 are derived from the an- 
gular lengths by assuming a distance of 315 pc and 500 pc to 
IC348 and OTC, respectively. The edge lengths for stars and 
BDs are listed individually for each cluster inside and outside 
the cluster full-width half maximum boundaries (also referred 
as cluster core radius R e ). In the IC348 cluster, the boundaries 
of 5'and 10' represent the cluster core radius R c and the limit 
at which the radial profi le merges with the background. These 
values are adopted from IMuench et al.l d2002l) where the radial 
distribution of sources is modelled. Similarly, a boundary of 
R C =1'2 is assumed for OT C based on the King Model fits from 
iHillenbrand & Hartmannl d!998l) . In IC348, the € MS r values for 
stars within R c =5'and outside vary by a factor of 1.9, whereas 
for substellar objects (denoted as BDs) they show similar val- 
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Fig. 1. IC348: a) Plot showing the variation of the Q parameter with magnitude of the point sources considered for the MST analysis. 
The solid line represents the curve for Flamingos K band data complete to K=18 mag. The dotted and dashed curves are for two 
control fields data obtained from 2MASS. b) The minimum spanning trees of stars (sources brighter than K=15 mag), c) The 
minimal spanning tree for BDs (sources fainter than K=15 mag). The inner and outer dotted circles have radii of 5'and 10'in b) and 
c). 




Fig. 2. Same as Fig. 1 for the Orion Trapezium Cluster. In this case the sources brighter than K=14 mag are considered as stars. The 
dotted circle marks a radius of 1 '2 representing the Trapezium cluster radius obtained by a King model fit. 



ues. In the case of OTC the ^mst values of stars change by a 
factor of 2.12 and those of BDs by 1.4 within and outside the 
assumed cluster core radius of R e = V.2. 



The above results are free from contamination due to binarity 
(as we excluded tree branches smaller than 150 AU) but include 
small N clusters. A total of 19 sources were excluded in the OTC 
sample that imitated binary lengths. The comparisions between 
inside and outside the cluster radius R e show that the stars are 
centrally concentrated in both IC348 and OTC clusters whereas 
the BDs are distributed uniformly. In IC348, the BDs are clearly 
homogeneo usly distributed within the radius of 10' reproducing 
the results of lMuench et all d2003l) . The MST variation of BDs in 
OTC may indicate weak concentration at the center, although the 
region examined here encompasses only the Trapezium cluster 
which is the central part of the much larger Orion Nebula Cluster 
(lHillenbrandl [l997). Therefore, these variations show, that un- 
like stars, the BDs are not limited by the cluster potential in both 



Table 1. Minimum Spanning Tree parameters for IC348 and 
OTC 





^MST 




f . 

L min 


/ 

1 max 


n, 


Survey 


10 3 AU 


10 3 AU 


10 3 AU 


10 3 AU 




IC348 


Stars (<5') 


7.8 


4.1 


1.4 


22.8 


135 


Stars (>5') 


14.7 


7.7 


0.8 


43.8 


238 


BDs (<5') 


9.6 


5.0 


1.1 


25.3 


98 


BDs (>5') 


9.5 


5.1 


1.0 


36.3 


536 


Orion Trapezium 


Stars (<1.2') 


2.5 


1.3 


0.4 


6.5 


228 


Stars (> 1.2') 


5.4 


3.4 


0.16 


26.1 


570 


BDs (<1.2') 


5.8 


3.4 


0.5 


13.7 


49 


BDs (>1.2') 


8.4 


5.0 


0.16 


26.9 


258 



IC348 and OTC. The substellar objects are rather uniformly dis- 
tributed in space imitating a fractal configuration. 
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4. Discussion 

The results of the analysis presented in the previous section has 
two significant parts: 

a) The BDs are spread uniformly in space and not limited 
by the cluster potential well whereas stars are bound by a well 
described cluster potential. This can be explained if BDs in the 
cluster are seen as an N particle system with higher initial ve- 
locities than the stars in random directions produced by mul- 
tiple interactions at early times. The net result will be a rela- 
tively uniform distribution of space positions and velocities in- 
stead of being constrained by the gradient of the cluster potential 
well. Such a scenario of multiple objects and chaotic interactions 
is well de s cribed by the numerical simulations such as that of 
iBate et al.1 d2003l) . However, this process is relevant during the 
early evolution of the cluster (within 10 5 yrs) thus setting up the 
"initial" conditions for BDs. Higher initial velocities are nec- 
essary to overcome the cluster potential gradient and produce 
the homogeneous distribution of spatial locations. The natural 
effect of lower mass objects acquirin g higher recoil velocities 
in comparision to higher mass objects ( Reip urth & Cla rke 2001; 
iRroupa & Bouvieri2 003) can result in BDs having the necessary 
ej ection velocities. This is i n accordance with the predictions 
of Kroupa ^fe Bouvier! d2003l) who compute higher dispersion in 
ejection velocities of 2-3 km s _1 for BDs. 

b) The BDs appear slightly more centrally concentrated in 
the 1 Myr old OTC than in the 3 Myr old IC348 cluster. While 
the BDs show a homogeneous distribution in space, they are still 
limited to the physical boundaries of the cluster or only slightly 
spread out. This is evident from the significant number of BDs in 
comparison to stars in both clusters (see last column of Table 1). 
It implies that the BDs have not traversed in space significantly 
away from the regions where they were born. The more uniform 
distribution in the 3 Myr old IC348 cluster with respect to the 
1 Myr old OTC imply that, with time, the BDs may indeed be 
moving away from the cluster center, but the net displacement 
is not significant. Assuming higher initial velocities for BDs, 
the cluster potential that operates on the larger scale, will re- 
duce their ove r all d isplacement away from the cluster center. 
iGoodwin et ail d2005l) examine the spatial distribution of BDs 
produced by the decay of small N stellar systems as expected 
from the embryo ejection scenario and show that significant spa- 
tial spread of BDs with respect to stars may occur only after 5 
Myr. Their predictions are consistent with the results of IC348 
and OTC. 

The above results are therefore cl ear signatures of the sce- 
nario of "ejected stellar embryos" (Reipurth & Clarkd 120011: 
iKrou pa & Bouvier 2003) as an important mechanism for the for- 
mation of brown dwarfs and free-floating planetary mass objects. 

Several studies have been made to investigate the differ- 
ences in spatial distribution of stars and BDs in star forming 
regions. In the Chamaeleon I cloud, the di stribution of BDs 
and lo w mass stars was found to be similar (lLopez Martf et all 
2004) and in the Lupus cloud BDs were found close to their 
birth sites dEd pez Martf et al J 120051) . In the Taurus molecular 
cloud, iBriceno et al.l d2002l) examined the distributions of near- 
est neighbours to objects with spectral indices <M6 and >M6 
and concluded that there were no significant differences in the 
spatial distribut ion of stars and BD s. Examining a smaller group 
of T-Tauri stars iMartfn et al.l d200ll) hinted at an anticorrel ation 
between density of stars and BDs in Taurus. ILuhmanl d2006l) ex- 
amined an area of 225 deg 2 covering almost the entire Taurus 
Molecular Cloud and found no significant differences in the spa- 
tial distributions of BDs and stars. 



The above studies were made because according to the 
"ejected embryos" scenario, BDs may be expected to be 
significantly displaced from their birth locations. For exam- 
ple,assuming a ID dispersion velocity of 1 kms 'they should 
have moved ~1 pc away (corresponding to 7'«6xR e ) from the 
cluster in lMyr old OTC or about 3pc away (corresponding to 
35'*7xR c ) in the IC348 cluster( if they do not experience any 
retarding force from the cluster potential). However, as evident 
from the data presented here, the BDs are indeed located within a 
region ~ 2xR e in both the clusters. Whether there is a significant 
concentration of BDs much farther away from the cluster cores 
in IC348 and OTC is an issue that needs investigation and is be- 
yond the scope of the data a nalysed here. H owever, following 
the example of Taurus, where ILuhmanl d2006l) examined the en- 
tire molecular cloud, significant populations of BDs away from 
the dense cores may not be expected. 

The major difference between previous studies and this work 
is the advantage of large number statistics and space densities of 
point sources found in IC348 and OTC. This number statistics 
has made it possible to effectively apply the MST technique to 
analyse the spatial distribution. In the absence of the MST anal- 
ysis, IC348 and OTC would imitate the resul ts of Taurus, Ch al 
or Lupus described above. Indeed in Fig. 2 of lLuhmanl (|2006) it 
can be seen that the BDs are sparsely populated while stars are 
relatively clustered. MST analysis for Taurus may also show the 
kind of variations found here, but the IC348 and OTC are clearly 
better data sets for effective application of the MST method 
which is "statistical" in nature. In conclusion, we have shown 
that the near-infrared observations of the IC348 and OTC em- 
bedded clusters are in accordance with the "ejected stellar em- 
bryos" scenario for the formation of substellar objects. Higher 
ejection velocities for substellar objects are essential but their 
spatial displacements from the birth sites may not be necessary 
to explain the observed data. 
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